Predicted genetic burden and frequency of phenotype-associated variants in the horse.

Disease-causing variants have been identified for less than 20% of suspected equine genetic diseases. Whole genome sequencing (WGS) allows rapid identification of rare disease causal variants. However, interpreting the clinical variant consequence is confounded by the number of predicted deleterious variants that healthy individuals carry (predicted genetic burden). Estimation of the predicted genetic burden and baseline frequencies of known deleterious or phenotype associated variants within and across the major horse breeds have not been performed. We used WGS of 605 horses across 48 breeds to identify 32,818,945 variants, demonstrate a high predicted genetic burden (median 730 variants/horse, interquartile range: 613-829), show breed differences in predicted genetic burden across 12 target breeds, and estimate the high frequencies of some previously reported disease variants. This large-scale variant catalog for a major and highly athletic domestic animal species will enhance its ability to serve as a model for human phenotypes and improves our ability to discover the bases for important equine phenotypes.

Genetic Variation and the Distribution of Variant Types in the Horse.

Genetic variation is a key contributor to health and disease. Understanding the link between an individual's genotype and the corresponding phenotype is a major goal of medical genetics. Whole genome sequencing (WGS) within and across populations enables highly efficient variant discovery and elucidation of the molecular nature of virtually all genetic variation. Here, we report the largest catalog of genetic variation for the horse, a species of importance as a model for human athletic and performance related traits, using WGS of 534 horses. We show the extent of agreement between two commonly used variant callers. In data from ten target breeds that represent major breed clusters in the domestic horse, we demonstrate the distribution of variants, their allele frequencies across breeds, and identify variants that are unique to a single breed. We investigate variants with no homozygotes that may be potential embryonic lethal variants, as well as variants present in all individuals that likely represent regions of the genome with errors, poor annotation or where the reference genome carries a variant. Finally, we show regions of the genome that have higher or lower levels of genetic variation compared to the genome average. This catalog can be used for variant prioritization for important equine diseases and traits, and to provide key information about regions of the genome where the assembly and/or annotation need to be improved.

Assessment of the FAM174A 11G allele as a risk allele for equine metabolic syndrome.

An 11G nucleotide repeat in the 3' UTR of FAM174A was recently postulated as a risk allele with a dominant mode of inheritance for equine metabolic syndrome (EMS) and laminitis status in Arabian horses. The objective of this project was to evaluate this hypothesis in a large and diverse across-breed population. A total of 301 ponies, 292 Morgans, 64 Arabians, 49 Tennessee Walking Horses and 59 Quarter Horses were genotyped for six observed G repeat alleles in the FAM174A 3' UTR. Phenotype data included laminitis status, baseline insulin, glucose, non-esterified fatty acids, triglycerides, adiponectin, leptin, ACTH, insulin and glucose post oral sugar test, and two proxies for insulin resistance. The 11G allele frequencies were 18.8, 6.9, 1.8, 0.2 and 0.0% in the Arabians, Tennessee Walkers, ponies, Morgans and Quarter Horses respectively. Association analyses between FAM174A genotype and EMS phenotypes, and between allele count and EMS phenotypes, identified no statistically significant associations. When a dominant effect for the 11G allele was evaluated, a statistically significant association with adiponectin levels was identified in the ponies, and pairwise comparisons revealed that the estimated marginal means were higher in ponies with the 11G allele vs. alternative alleles (i.e. the allele had a protective effect). In conclusion, our data do not support the FAM174A 11G allele as a risk allele for EMS in our studied breeds.

Associations between endocrine disrupting chemicals and equine metabolic syndrome phenotypes.

Equine Metabolic Syndrome (EMS) is characterized by abnormalities in insulin regulation, increased adiposity and laminitis, and has several similarities to human metabolic syndrome. A large amount of environmental variability in the EMS phenotype is not explained by commonly measured factors (diet, exercise, and season), suggesting that other environmental factors play a role in EMS development. Endocrine disrupting chemicals (EDCs) are associated with metabolic syndrome and other endocrine abnormalities in humans. This led us to hypothesize that EDCs are detectable in horse plasma and play a role in the pathophysiology of EMS. EDCs acting through the aryl hydrocarbon and estrogen receptors, were measured in plasma of 301 horses from 32 farms. The median (range) TEQ (2,3,7,8-TCDD equivalent) and EEQ (17ß-estradiol equivalent) were 19.29 pg/g (0.59-536.36) and 10.50 pg/ml (4.35-15000.00), respectively. TEQ was negatively associated with plasma fat extracted and batch analyzed. EEQ was positively associated with pregnancy and batch analyzed, and negatively associated with being male and superfund score ≤100 miles of the farm. Of particular interest, serum glucose and insulin, glucose and insulin post oral sugar challenge, and leptin concentrations were associated with EEQ, and serum triglyceride concentration was associated with TEQ. Overall, we demonstrated that EDCs are present in the plasma of horses and may explain some of the environmental variability in measured EMS phenotypes. This is the first example of EDCs being associated with clinical disease phenotype components in domestic animals.

Heritability of metabolic traits associated with equine metabolic syndrome in Welsh ponies and Morgan horses.

BACKGROUND: Equine metabolic syndrome (EMS) is a complex clinical disorder with both environmental and genetic factors contributing to EMS phenotypes. Estimates of heritability determine the proportion of variation in a trait that is attributable to genetics. OBJECTIVES: To provide heritability estimates for nine metabolic traits associated with EMS in two high-risk breeds. STUDY DESIGN: Retrospective cohort study. METHODS: High-density single-nucleotide polymorphism (SNP) genotype data was used to estimate the heritability (h2 SNP ) of nine metabolic traits relevant to EMS in a cohort of 264 Welsh ponies and 286 Morgan horses. Traits included measurements of insulin, glucose, non-esterified fatty acids (NEFA), triglycerides, leptin, adiponectin, ACTH, and glucose (GLU-OST) and insulin (INS-OST) following an oral sugar challenge. RESULTS: In Welsh ponies, seven of the nine traits had statistically significant h2 SNP estimates that were considered moderately to highly heritable (h2 SNP >0.20) including: triglycerides (0.313; s.e. = 0.146), glucose (0.408; s.e. = 0.135), NEFA (0.434; s.e. = 0.136), INS-OST (0.440; s.e. = 0.148), adiponectin (0.488; s.e. = 0.143), leptin (0.554; s.e. = 0.132) and insulin (0.808; s.e. = 0.108). In Morgans, six of the nine traits had statistically significant h2 SNP estimates that were also determined to be moderately to highly heritable including: INS-OST (0.359; s.e. = 0.185), leptin (0.486; s.e. = 0.177), GLU-OST (0.566 s.e. = 0.175), insulin (0.592; s.e. = 0.195), NEFA (0.684; s.e. = 0.164), and adiponectin (0.913; s.e. = 0.181). MAIN LIMITATIONS: Insufficient population size may have limited power to obtain statistically significant h2 SNP estimates for ACTH (both breeds), glucose and triglycerides in Morgans and GLU-OST in Welsh ponies. CONCLUSIONS: This study provides the first concrete evidence of a genetic contribution to key phenotypes associated with EMS. Eight of these nine traits had moderate to high h2 SNP estimates in this cohort. These data demonstrate that continued research for identification of the genetic risk factors for EMS phenotypes within and across breeds is warranted.

SNP-based heritability and genetic architecture of tarsal osteochondrosis in North American Standardbred horses.

Osteochondrosis is a common developmental orthopedic disease characterized by a failure of endochondral ossification. Standardbred horses are recognized as being predisposed to tarsal osteochondrosis. Prior heritability estimates for tarsal osteochondrosis in European Standardbreds and related trotting breeds have been based on pedigree data and range from 17-29%. Here, we report on genetic architecture and heritability based on high-density genotyping data in a cohort of North American Standardbreds (n = 479) stringently phenotyped for tarsal osteochondrosis. Whole-genome array genotyping data were imputed to ~2 million single nucleotide polymorphisms (SNPs). SNP-based heritability of osteochondrosis in this population was explained by 2326 SNPs. The majority of these SNPs (86.6%) had small effects, whereas fewer SNPs had moderate or large effects (10% and 2.9% respectively), which is consistent with a polygenic/complex disease. Heritability was estimated at 0.24 ± 0.16 using two methods of restricted maximum likelihood analysis, as implemented in gcta (with and without a weighted relatedness matrix) and ldak software. Estimates were validated using bootstrapping. Heritability estimates were within the range previously reported and suggest that osteochondrosis is moderately heritable but that a significant portion of disease risk is due to environmental factors and/or genotype × environment interactions. Future identification of the genes/variants that have the most impact on disease risk may allow early recognition of high-risk individuals.

Canine NAPEPLD-associated models of human myelin disorders.

Canine leukoencephalomyelopathy (LEMP) is a juvenile-onset neurodegenerative disorder of the CNS white matter currently described in Rottweiler and Leonberger dogs. Genome-wide association study (GWAS) allowed us to map LEMP in a Leonberger cohort to dog chromosome 18. Subsequent whole genome re-sequencing of a Leonberger case enabled the identification of a single private homozygous non-synonymous missense variant located in the highly conserved metallo-beta-lactamase domain of the N-acyl phosphatidylethanolamine phospholipase D (NAPEPLD) gene, encoding an enzyme of the endocannabinoid system. We then sequenced this gene in LEMP-affected Rottweilers and identified a different frameshift variant, which is predicted to replace the C-terminal metallo-beta-lactamase domain of the wild type protein. Haplotype analysis of SNP array genotypes revealed that the frameshift variant was present in diverse haplotypes in Rottweilers, and also in Great Danes, indicating an old origin of this second NAPEPLD variant. The identification of different NAPEPLD variants in dog breeds affected by leukoencephalopathies with heterogeneous pathological features, implicates the NAPEPLD enzyme as important in myelin homeostasis, and suggests a novel candidate gene for myelination disorders in people.

A highly prevalent equine glycogen storage disease is explained by constitutive activation of a mutant glycogen synthase.

BACKGROUND: Equine type 1 polysaccharide storage myopathy (PSSM1) is associated with a missense mutation (R309H) in the glycogen synthase (GYS1) gene, enhanced glycogen synthase (GS) activity and excessive glycogen and amylopectate inclusions in muscle. METHODS: Equine muscle biochemical and recombinant enzyme kinetic assays in vitro and homology modelling in silico, were used to investigate the hypothesis that higher GS activity in affected horse muscle is caused by higher GS expression, dysregulation, or constitutive activation via a conformational change. RESULTS: PSSM1-affected horse muscle had significantly higher glycogen content than control horse muscle despite no difference in GS expression. GS activity was significantly higher in muscle from homozygous mutants than from heterozygote and control horses, in the absence and presence of the allosteric regulator, glucose 6 phosphate (G6P). Muscle from homozygous mutant horses also had significantly increased GS phosphorylation at sites 2+2a and significantly higher AMPKα1 (an upstream kinase) expression than controls, likely reflecting a physiological attempt to reduce GS enzyme activity. Recombinant mutant GS was highly active with a considerably lower Km for UDP-glucose, in the presence and absence of G6P, when compared to wild type GS, and despite its phosphorylation. CONCLUSIONS: Elevated activity of the mutant enzyme is associated with ineffective regulation via phosphorylation rendering it constitutively active. Modelling suggested that the mutation disrupts a salt bridge that normally stabilises the basal state, shifting the equilibrium to the enzyme's active state. GENERAL SIGNIFICANCE: This study explains the gain of function pathogenesis in this highly prevalent polyglucosan myopathy.

Genetic risk factors for insidious equine recurrent uveitis in Appaloosa horses.

Appaloosa horses are predisposed to equine recurrent uveitis (ERU), an immune-mediated disease characterized by recurring inflammation of the uveal tract in the eye, which is the leading cause of blindness in horses. Nine genetic markers from the ECA1 region responsible for the spotted coat color of Appaloosa horses, and 13 microsatellites spanning the equine major histocompatibility complex (ELA) on ECA20, were evaluated for association with ERU in a group of 53 Appaloosa ERU cases and 43 healthy Appaloosa controls. Three markers were significantly associated (corrected P-value <0.05): a SNP within intron 11 of the TRPM1 gene on ECA1, an ELA class I microsatellite located near the boundary of the ELA class III and class II regions and an ELA class II microsatellite located in intron 1 of the DRA gene. Association between these three genetic markers and the ERU phenotype was confirmed in a second population of 24 insidious ERU Appaloosa cases and 16 Appaloosa controls. The relative odds of being an ERU case for each allele of these three markers were estimated by fitting a logistic mixed model with each of the associated markers independently and with all three markers simultaneously. The risk model using these markers classified ~80% of ERU cases and 75% of controls in the second population as moderate or high risk, and low risk respectively. Future studies to refine the associations at ECA1 and ELA loci and identify functional variants could uncover alleles conferring susceptibility to ERU in Appaloosa horses.

Coat color genotypes and risk and severity of melanoma in gray quarter horses.

BACKGROUND: Both graying and melanoma formation in horses have recently been linked to a duplication in the STX17 gene. This duplication, as well as a mutation in the ASIP gene that increases MC1R pathway signaling, affects melanoma risk and severity in gray horses. OBJECTIVE: To determine if melanoma susceptibility in gray Quarter Horses (QH) is lower than gray horses from other breeds because of decreased MC1R signaling resulting from a high incidence of the MC1R chestnut coat color allele in the QH population. ANIMALS: A total of 335 gray QH with and without dermal melanomas. METHODS: Blood or hair root samples were collected from all horses for DNA extraction and genotyping for STX17, ASIP, and MC1R genotypes. Age, sex, and external melanoma presence and grade were recorded. The effect of age and genotype on melanoma presence and severity was evaluated by candidate gene association. RESULTS: Melanoma prevalence (16%) and grade (0.35) in this QH cohort was lower than that reported in other breeds. Age was significantly associated with melanoma prevalence (P = 5.28 × 10(-11)) and severity (P = 2.2 × 10(-13)). No significant effect of MC1R genotype on melanoma prevalence or severity was identified. An effect of ASIP genotype on melanoma risk was not detected. Low STX17 homozygosity precluded evaluation of the gray allele effect. CONCLUSION AND CLINICAL IMPORTANCE: Melanoma prevalence and severity is lower in this population of gray QH than what is reported in other breeds. This could be because of the infrequent STX17 homozygosity, a mitigating effect of the MC1R mutation on ASIP potentiation of melanoma, other genes in the MC1R signaling pathway, or differences in breed genetic background.

Equine developmental orthopaedic diseases--a genome-wide association study of first phalanx plantar osteochondral fragments in Standardbred trotters.

Palmar/plantar osteochondral fragments (POF) in fetlock joints commonly affect and influence the athletic performance of horses. In this study, we used the Equine SNP50 BeadChip® to perform a genome-wide association study of metatarsophalangeal POF in 176 Norwegian Standardbred trotter yearlings. Putative quantitative trait loci (QTL) for medial and/or lateral POF, and medial POF only were identified on ECA1, 2, 7, 9 and 31, whereas for lateral POF, only on ECA7, 11, 27 and X. The moderate number of QTL evidences a complex inheritance and suggests various genes controlling POF development in medial and lateral locations.

Genetic mapping of recurrent exertional rhabdomyolysis in a population of North American Thoroughbreds.

Recurrent exertional rhabdomyolysis is a heritable disorder that results in painful skeletal muscle cramping with exercise in up to 10% of all Thoroughbred racehorses. Here, we report a genome-wide association study with 48 282 SNPs analyzed among 48 case and 37 control Thoroughbreds. The most significant SNPs spanned approximately 13 Mb on ECA16, and the P-value of the most significant SNP after correcting for population structure was 8.0 × 10(-6) . This region on ECA16 was further evaluated by genotyping 247 SNPs in both the initial population and a second population of 34 case and 98 control Thoroughbreds. Several SNPs across the 13-Mb region on ECA16 showed significance when each population was analyzed separately; however, the exact positions of the most significant SNPs within this region on ECA16 varied between populations. This variability in location may be attributed to lack of power owing to insufficient sample sizes within each population individually, or to the relative distribution of long, conserved haplotypes, characteristic of the Thoroughbred breed. Future genome-wide association studies with additional horses would likely improve the power to resolve casual loci located on ECA16 and increase the likelihood of detecting any additional loci on other chromosomes contributing to disease susceptibility.

Epidemiological and genetic study of exertional rhabdomyolysis in a Warmblood horse family in Switzerland.

REASONS FOR PERFORMING STUDY: Exertional rhabdomyolysis (ER) and its familial basis in Warmblood horses is incompletely understood. OBJECTIVES: To describe the case details, clinical signs and management of ER-affected Warmblood horses from a family with a high prevalence of ER, to determine if histopathological signs of polysaccharide storage myopathy (PSSM) and the glycogen synthase (GYS1) mutation are associated with ER in this family, and to investigate potential risk factors for development of ER. METHODS: A family consisting of a sire with ER and 71 of his descendants was investigated. History of episodes of ER, husbandry, feeding and use was assessed by interviewing horse owners using a standardised questionnaire. All horses were genotyped for GYS1. In 10 ER-affected horses, muscle histopathology was evaluated. RESULTS: Signs of ER were reported in 39% of horses and 51% of the entire family possessed the GYS1 mutation. Horses possessing the GYS1 mutation had a 7.1-times increased risk for developing ER compared to those with the normal genotype (95% confidence interval [CI] 2.37-21.23, P = 0.0005). All muscle samples from horses in the family with ER showed polysaccharide accumulation typical for PSSM, amylase-resistant in 9/10 cases. There was evidence (odds ratio 5.6, CI 1.00-31.32, P = 0.05) that fat or oil feeding improved clinical signs of ER. No other effects of environmental factors associated with clinical signs of ER were identified. CONCLUSION AND POTENTIAL RELEVANCE: PSSM associated with the GYS1 mutation is one identifiable cause of ER in Warmblood horses. Signs of ER are not always manifest in GYS1 positive horses and there are also other causes for ER in Warmblood horses. Breeding animals with the GYS1 mutation results in a high prevalence of ER due to its dominant mode of inheritance.

Genome-wide association analysis of osteochondrosis of the tibiotarsal joint in Norwegian Standardbred trotters.

Osteochondrosis (OC), a disturbance in the process of endochondral ossification, is by far the most important equine developmental orthopaedic disease and is also common in other domestic animals and humans. The purpose of this study was to identify quantitative trait loci (QTL) associated with osteochondrosis dissecans (OCD) at the intermediate ridge of the distal tibia in Norwegian Standardbred (SB) using the Illumina Equine SNP50 BeadChip whole-genome single-nucleotide polymorphism (SNP) assay. Radiographic data and blood samples were obtained from 464 SB yearlings. Based on the radiographic examination, 162 horses were selected for genotyping; 80 of these were cases with an OCD at the intermediate ridge of the distal tibia, and 82 were controls without any developmental lesions in the joints examined. Genotyped horses descended from 22 sires, and the number of horses in each half-sib group ranged from 3 to 14. The population structure necessitated statistical correction for stratification. When conducting a case-control genome-wide association study (GWAS), mixed-model analyses displayed regions on chromosomes (Equus callabus chromosome - ECA) 5, 10, 27 and 28 that showed moderate evidence of association (P ≤ 5 × 10(-5); this P-value is uncorrected i.e. not adjusted for multiple comparisons) with OCD in the tibiotarsal joint. Two SNPs on ECA10 represent the most significant hits (uncorrected P=1.19 × 10(-5) in the mixed-model). In the basic association (chi-square) test, these SNPs achieved statistical significance with the Bonferroni correction (P=0.038) and were close in the permuted logistic regression test (P=0.054). Putative QTL on ECA 5, 10, 27 and 28 represent interesting areas for future research, validation studies and fine mapping of candidate regions. Results presented here represent the first GWAS of OC in horses using the recently released Illumina Equine SNP50 BeadChip.

Estimated prevalence of the Type 1 Polysaccharide Storage Myopathy mutation in selected North American and European breeds.

The GYS1 gene mutation that is causative of Type 1 Polysaccharide Storage Myopathy (PSSM) has been identified in more than 20 breeds of horses. However, the GYS1 mutation frequency or Type 1 PSSM prevalence within any given breed is unknown. The purpose of this study was to determine the frequency of the GYS1 mutation and prevalence of genetic susceptibility to Type 1 PSSM in selected breeds from Europe and North America. The GYS1 mutation was detected in 11 breeds, including, in order of increasing allele frequency, Shires, Morgans, Appaloosas, Quarter Horses, Paints, Exmoor Ponies, Saxon-Thuringian Coldbloods, South German Coldbloods, Belgians, Rhenish German Coldbloods and Percherons. The prevalence of genetic susceptibility to Type 1 PSSM in these breeds varied from 0.5% to 62.4%. The GYS1 mutation was not found in the sampled Thoroughbreds, Akhal-Tekes, Connemaras, Clydesdales, Norwegian Fjords, Welsh Ponies, Icelandics, Schleswig Coldbloods or Hanoverians, but failure to detect the mutation does not guarantee its absence. This knowledge will help breed associations determine whether they should screen for the GYS1 mutation and will alert veterinarians to a possible differential diagnosis for muscle pain, rhabdomyolysis or gait abnormalities.

Presence of the glycogen synthase 1 (GYS1) mutation causing type 1 polysaccharide storage myopathy in continental European draught horse breeds.

The purpose of this study was to determine which continental European draught horse breeds harbour a mutation in the glycogen synthase 1 gene (GYS1) that is known to be responsible for type 1 polysaccharide storage myopathy in quarter horses and North American draught horses. Of a non-random selection of continental European draught horses belonging to 13 breeds, 62 per cent (250 of 403) tested were found to carry the mutant allele. The horses were located in Belgium, France, Germany, The Netherlands, Spain and Sweden. The mutation was identified in animals from each of the breeds examined. In the breeds in which more than 15 animals were available for testing, the highest percentages of GYS1-positive horses were found in the Belgian trekpaard (92 per cent; 35 of 38 horses tested), Comtois (80 per cent; 70 of 88), Netherlands trekpaard (74 per cent; 17 of 23), Rheinisch-Deutsches kaltblut (68 per cent; 30 of 44) and Breton (64 per cent; 32 of 51).

Genome sequence, comparative analysis, and population genetics of the domestic horse.

We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

A glycogen synthase 1 mutation associated with equine polysaccharide storage myopathy and exertional rhabdomyolysis occurs in a variety of UK breeds.

REASONS FOR PERFORMING STUDY: A glycogen synthase (GYS1) mutation has been described in horses with histopathological evidence of polysaccharide storage myopathy (PSSM) in the USA. It is unknown whether the same mutation is present in horses from the UK. OBJECTIVES: To determine whether the GYS1 mutation occurs in UK horses with histopathological evidence of PSSM and exertional rhabdomyolysis. HYPOTHESIS: The R309H GYS1 mutation is present in a variety of UK horse breeds and that the mutation is commonly associated with exertional rhabdomyolysis. METHODS: DNA was extracted from 47 muscle or blood samples from UK horses with histories of exertional rhabdomyolysis in which muscle biopsy diagnosis had been pursued. The proportions of GYS1 mutation positive cases were compared among histopathologically defined groups. In addition, breeds that carried the GYS1 mutation were identified from a total of 37 grade 2 (amylase-resistant) PSSM cases. RESULTS: Of 47 horses with exertional rhabdomyolysis in which a muscle biopsy diagnosis was pursued, 10 (21%) carried the GYS1 mutation. The mutation was only found in horses with grade 2 PSSM (i.e. not in horses with normal, idiopathic myopathy or grade 1 PSSM biopsy samples). In total, the GYS1 mutation was found in 24/37 (65%) of grade 2 PSSM cases. A variety of breeds, including Quarter Horse, Appaloosa, Warmblood, Connemara-cross, Cob, Polo Pony and Thoroughbred cross carried the mutation. CONCLUSIONS: The GYS1 mutation is an important cause of exertional rhabdomyolysis of UK horse breeds but does not account for all forms of PSSM. POTENTIAL RELEVANCE: Genotyping is recommended in cases of exertional rhabdomyolysis, prior to or in combination with, muscle biopsy. However a significant proportion of horses with histopathological evidence of PSSM and/or exertional rhabdomyolysis have different diseases.

Comparative skeletal muscle histopathologic and ultrastructural features in two forms of polysaccharide storage myopathy in horses.

Polysaccharide storage myopathy (PSSM) has been found in more than 35 different horse breeds through identification of abnormal storage of polysaccharide in muscle biopsies. A dominant mutation in the glycogen synthase 1 gene (GYS1) accounts for a substantial proportion of PSSM cases in at least 17 breeds, including Quarter Horses, but some horses diagnosed with PSSM by muscle histopathologic analysis are negative for the mutation. We hypothesized that a second distinct form of glycogen storage disease exists in GYS1-negative horses with PSSM. The objectives of this study were to compare the histopathologic features, ultrastructure of polysaccharide, signalment, history, and presenting complaints of GYS1-negative Quarter Horses and related breeds with PSSM to those of GYS1-positive horses with PSSM. The total histopathologic score in frozen sections of skeletal muscle stained with hematoxylin and eosin, periodic acid Schiff (PAS) and amylase-PAS stains from 53 GYS1-negative horses did not differ from that of 52 GYS1-positive horses. Abnormal polysaccharide was fine granular or homogenous in appearance (49/53; 92%), often amylase-sensitive (28/53; 53%), more commonly located under the sarcolemma, and consisting of beta glycogen particles in GYS1-negative horses. However, in GYS1-positive horses, abnormal polysaccharide was usually coarse granular (50/52; 96%), amylase-resistant (51/52; 98%), more commonly cytoplasmic, and consisting of beta glycogen particles or, in some myofibers, filamentous material surrounded by beta glycogen particles. Retrospective analysis found that GYS1-negative horses (n = 43) were younger at presentation (4.9 +/- 0.6 years vs. 6.7 +/- 0.3 years for GYS1-positive horses) and were more likely to be intact males than GYS1-positive horses (n = 160). We concluded that 2 forms of PSSM exist and often have distinctive abnormal polysaccharide. However, because evaluation of the histologic appearance of polysaccharide can be subjective and affected by age, the gold standard for diagnosis of PSSM at present would appear to be testing for the GYS1 mutation followed by evaluating muscle biopsy for characteristic abnormal polysaccharide in those horses that are negative for the mutation.